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ABSTRACT 

The Gemini Planet Imager (GPI) is a high performance adaptive optics system being designed 
and built for the Gemini Observatory. GPI is optimized for high contrast imaging, combining precise 
and accurate wavefront control, diffraction suppression, and a speckle-suppressing science camera with 
integral field and polarimetry capabilities. The primary science goal for GPI is the direct detection and 
characterization of young, Jovian-mass exoplanets. For plausible assumptions about the distribution 
of gas giant properties at large semi-major axes, GPI will be capable of detecting more than 10% 
of gas giants more massive than 0.5 Mj around stars younger than 100 Myr and nearer than 75 
parsecs. For systems younger than 1 Gyr, gas giants more massive than 8 Mj and with semi-major 
axes greater than 15 AU are detected with completeness greater than 50%. A survey targeting young 
stars in the solar neighborhood will help determine the formation mechanism of gas giant planets by 
studying them at ages where planet brightness depends upon formation mechanism. Such a survey 
will also be sensitive to planets at semi-major axes comparable to the gas giants in our own solar 
system. In the simple, a nd id ealized, situation in which plane ts formed by either the "hot-start" 
model of iBurrows et al.l ()2003l ) or the core accretion model of iMarley et al.l ([2007) , a few tens of 
detected planets are sufficient to distinguish how planets form. 
Subject headings: Extrasolar planets; high contrast imaging 



1. INTRODUCTION 

The past decade has seen extraordinary progress in 
detection of exoplanets with the deployment of cus- 
tom designed sys tems that use the Doppler method, 
such as HARPS (IMavor et al.l I2003D. and the transit 
method, such as Kepler (jBorucki et al.H2003h . Exoplan- 
etar y systems are also being discovered us ing microlens- 
ing (jGould et al.l |20H iSumi et al.l [20ToT ) . No planets 
are confirmed to have been detected using astrometry, 
though there are disputed clai ms of some detections, in- 
cluding a compan i on to VB 10 ([Pravdo fc Shaklanl 120091 : 
IBean et al.l 120101: lAnglada-Escude et al.l 120100 . Even 
so, the upcoming Gaia mission is expected to have 
astrometric precision of ~ 8/ias, and thus be able 
to astrometrically discover the hundreds of exoplanets 
and exoplanetary systems expected to have astromet- 
ric signals > 30 uas around stars within 200 parsecs 
(jCasertano et al.1 12008). Polarimetric detection of HP 
18973 3 b has also been claim ed by iBerdvugina et al.l 
(2008]), but lWiktorowic3 p009) was unable to reproduce 
the detection. 

These methods are indirect, and provide different in- 
formation about the planets that they detect. From ra- 
dial velocity surveys alone, the period, eccentricity and 
a minimum mass may be measured for detected planets. 
For planets that also transit their hosts, it is possible to 
learn the inclination of the orbit and planet radius, given 
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simple assumptions about the system and knowledge 
of the stellar mass-radius rela tion, which then give the 
mass and density of the planet (jSeager fc Mallen-Ornelas 
. 2003 ;) . Mid-IR exoplanetary light has been detected 
in se c ondary eclipses by Sp it zer dCharbonneau et all 
20051: iGrillmair et all 120071: iRichardson et al.l 120071: 



Fressin et al.l 120101 : lO'Donovan et alll2010T ). though in- 
formation about the planet spectrum itself is limited by 
t he photon sho t noise of the primary. 

iMarcv et all (|2008f) lists 228 extrasolar planets, with 
5% of targeted stars possessing massive planets, and 
shows that a diversity of exoplanet systems exists. 
Though the number of radial velocity confirmed plan- 
ets has nearly doubled since0 and many more planet 
candidates exist in t he first four months of Kepler data 
(|Borucki et al.ll201l|) . radial velocity and transit surveys 
leave several long-standing questions about planetary 
systems unanswered. How do planets form? Is the solar 
system typical? What is the abundance of solar-like sys- 
tems? These surveys also raise new questions, including 
what produces the dynamical diversity in exoplanetary 
systems? 

The prospects for indirect planet detection techniques 
alone to answer these questions is limited. With the ex- 
ception of microlensing, the effectiveness of these tech- 
niques at detecting exoplanets decreases at large semi- 
major axes. For a reliable detection with radial velocity, 
a significant fraction of an orbital period must elapse. As 
such, radial velocity surveys are only now reaching the 
precision and lifetime necessary to detect Jupiter, and 
thus do not yet constrain the frequency of solar system 
analogs. For example, i n the Keck radial y elocity search 
that began in 1996 July (B utler et al.ll2006Tl only planets 



5 The Exoplanet Orbit Dat abase at http://exo planets.org lists 
429 planets as of March 2010 HWright et al.H201(JI ). 
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with a < 6 AU have completed one orbit. The median 
semi-major axis of known exoplanet orbits is approxi- 
mately 1 AU, a nd HP 190984 b i s the only planet with 
with a > 6 AU (Santos et al. 2010). As the time baseline 
of surveys grows, the semi-major axes probed increase ac- 
cording to a — P 2 / 3 , meaning improving the statistics at 
large semi-major axes will be challenging. Doppler sur- 
veys have also been focused on a small range of stellar 
masses and ages, as precise radial velocity measurement 
requires strong stellar spectral features and low activ- 
ity. This has limited studies of exoplanet trends with 
stellar properties beyond basic correlations between host 
star mass and metallicity and the probability o f hosting 
a planet within a few AU (| Johnson et aLll2010D . 

Direct imaging constitutes the next step in character- 
izing other planetary systems. Many previous attempts 
to directly image substellar companions to stars have 
already placed stringent constraints on the presence 

of gia nt planets on wide or bits (ISchroeder et all 

200a IMasciadri et all [2005t iKasper et all 12007 



Biller et all l2007t iLafreniere et all l2007bl: lApai et al l 
2008tlHeinze et alJl^oToTlLeconte et al.ll2010l ). There are 



also p rojects already underwa y, such as NICI (|Liu et all 
I2010T ) and SEEDS (jTamural I2009D. or that are near 
deploying, such as SPHERE fBeuzit et all 120101 ) and 
GPL Thus far, planet detections through direct imaging 
have been limited by difficulty in achieving high contrast 
at small angular separations, but instruments are now 
reaching the threshold at which planet detection by 
direct imaging is promising. Direct imaging is highly 
complementary to Doppler surveys. Unlike Doppler 
searches, direct imaging is sensitive at large semi-major 
axes, as planets can be found without waiting for an 
orbit to complete — a condition that renders Doppler 
detection of planets on Neptune-like orbits with semi- 
major axes of 30 AU and periods of 160 yr impractical. 
Fourier decomposition, which underlies Doppler and 
astrometric detection, is also subject to aliasing and 
beat phenomena, and suffers confusion when multiple 
planets are present. A direct search that can probe 
beyond 5 AU will thus bring statistical significance to 
these studies. 

The goal of direct detection is to spatially separate 
the exoplanet light from that of its primary. This affords 
access to exoplanet atmospheres, which yields fundamen- 
tal information including effective temperature, gravity, 
atmospheric composition and abundances, orbital mo- 
tion, and perhaps even weather and planetary spin (via 
polarization a ssociated with rotation-induced oblateness; 
iMarlevI [20101 Direct imaging may also reveal the role 
of giant impacts, though 10 to 100 planets ~ 10 Myr old 
may need to be discovered in order to expect to find a 
single gas giant post impact (|Anic et al.l 120071 ) . or lead 
to discovery of time dependent atmospheric phenomena. 

Most importantly, by probing large semi-major axes, 
direct imaging will see gas giants beyond the "snow line," 
which is where they are expected to form, and extend 
out to the greatest distances at which giant planets can 
form. The location of the region of interest depends on 
at least two competing factors: time-scales for planet 
building and the availability of raw material. Dynam- 
ical and viscous time scales in the disk are shorter at 
small radii, while for typical surface-density laws the 



amount of mass increases with radius, with a jump in 
the abundance of solid material beyond the "snow line" 
where ices condense. This change in the surface den- 
sity of solid material occurs at 2.7 AU in the Hayashi 
model (lHayashi l 119811 ). The lo cation of this boundar y 
depends on the disk structure (jSasselov fc Lecarl 120001 ). 
but for solar type stars, the zone of interest is beyond 
that which is readily probed by the Doppler method. 
The discovery of giant planets far beyond the snow line 
would tend to favor theories of planet formation by grav- 
itational instability over solid core condensation and ac- 
cretion. At larger orbital radii (> 20-30 AU), gas- 
cooling times become shorter than the Keplerian shear- 
ing time (e.g., iKratter et al.ll20iol) — a n ecessary condi- 
tion for runaway gravitational instabi lity (Gammiel l2001l : 
Uohnson fc Gam mie 2003; Boss 2002j) — while solid core 
growth by collisional coag ulation of planetesima ls pro- 
ceeds prohibitively slowly ([Goldreich et al.l 120041 ) . Even 
so, new ideas about coupling between migration and core 
accretion suggest significantly increased growth, allow- 
ing the formation o f planets well beyond the snow line 
(jLevison et al.ll20Toh . The structure of our own solar sys- 
tem implies that a full picture of planet formation can- 
not be constructed without reaching out to 30 or 40 AU. 
Millimeter observations of T Tauri disks sup port this, as 
typic al disk radii fall in the range 50-100 AU (|Isella et all 
[20091) . 

The initial conditions, composition, and equation of 
state all influence evolution of young gas giant plan- 
ets. Gas giants will take tens to hundreds of mil- 
lions of years to "forget" their post-formation entropy, 
meaning temperatures and luminosities of you ng plan- 
ets will reveal their past ()Marlev et all I2007D . Stel- 
lar properties, in particular mass and metallicity, will 
also influence planet formation. The observed planet- 
metallicity correlation in Doppler planets has been in- 
terpreted as supporting core accretion as the forma- 
tion mechanism of gas giants ()Fischer fc ValentH 120051 : 
Uohnson et al.ll2~010D . Metal-rich disks are exp ected to 
have longer lifetimes (lErcolano fc C larke 201(J and en- 
hanced clump formation (| Johansen et al.l [20091 ) . all in- 
creasing the likelihood of planet formation for those 
formed by core accretion. Formation by disk instabil- 
ity is not e xpected to be sensitive to stellar metallicity 
(Boss 2002). This has led to the suggestion that low 
metallicity stars will preferentially form planets by disk 
instability on wide orbits, while high metallicity stars will 
host planets on short period orbits formed by co re accre- 
tion (|Bolevl 120091: iDodson-Robinson et alll2009D . Stellar 
metallicity is also expected to influence planet appear- 
ance. Planets formed via core accretion ar e predicted to 
them selves have super-st ellar metallicity ([Pollack et all 
[19961 : iFortnev et allboOSl ). Though gas giants more mas- 
sive than ~ 5Mj formed by disk instability may also 
have super-stellar atmospheric metallicities, less mas- 
sive gas giants should ha ve metal depleted atmospheres 
(|Helled fc S chubert 2008). Exoplanet infrared colors are 
predicted to be sensitive to atmospheric metal abun- 
dances, so direct imaging and spectrosco py will reveal 
inform ation about a planet's formation (jFortnev et all 
[2001) . 

Another reason to image the outer regions of extra- 
solar systems is to probe them for vestiges of planetary 
migration. Ninety percent of the Doppler sample con- 
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sists of massive planets with a < 3 AU, suggesting that 
they migrated inwards to their present locations. A vari- 
ety of mechanisms may drive orbital evolution; the tidal 
gravitatio nal interaction between the planet and a vis- 
cous disk ( Goldreic h fe TremaindlT980f ) , the gravitational 
inte raction between tw o or more Jupiter mass plan- 
ets (jRasio fe Ford [1996:) . and t he interaction betwe en a 
planet and a planetesimal disk ()Murrav et al.lfl998h . It 
is energetically favorable for a Keplerian disk to evolve 
by trans porting mass inward and ang ular momentum 
outward (|Lvnden-Bell fc Pringld 119741 ) . Consequently, 
inward planetary drift appears ine vitable, and this is 
what is found in certain simulations ([Trilling et al1l2002t 
lArmitage et al.ll200a iMatsuvama et al.ll2003[ ). However, 
if planets form while the disk is being dispersed, or if 
multiple planets are present, outward migration can also 
occur. In a system consisting initially of two Jupiter- 
like planets, a dynamical instability may eject one planet 
while the other is left in a tight, eccentric orbit. The sec- 
ond planet is not always lost; the observed Doppler exo- 
planct eccentricity distribution can be reproduced if the 
51 Pegasi systems are formed by planet-planet scattering 
events and the second planet typical ly remains bound in 
a wide (a > 20 AU), eccentric orbit dRasio fc FordllT996l : 
iMarzari fc Weidenschillinl[2002l: IVeras et al.ll2009l) . Di- 
vergent migration of pairs of Jupiter-mass planets within 
viscous disks leads to mutual resonance crossings and ex- 
citation of orbital eccentricities such that the resultant 
ellipticities are inve rsely correlated with planet masses 
([Chiang et al.l 120021 ) . Given decreasing disk viscosity 
with radius and the consequent reduction in planetary 
mobility with radius, we expect eccentricities to decrease 
with radius, perhaps s harply if the ma gneto-rotational 
instability is invoked ([Sano et al.l [20001 . By contrast, 
excitation of eccentricity by disk-planet interactions re- 
quires no additional planet to explai n the ellipticities 
of cu rrently known solitary planets (|Goldreich fc Saril 
2003). Clearly, observations of the incidence, mass, 
and eccentricity distributions of multiple planet systems 
would sharpen ideas regarding how planetary orbits are 
sculpted. 

There have already been a number of planets and 
planet candidates discovered by direct imaging, in clud- 
ing beta Pic b (8 AU/8 Mj; lLagrange et ail 120091) . the 
upper Sco object 1RX S J160929.1-210524 (330 AU/8 Mj 
iLafreniere etal.1 120081 ). HR 8799 b, c, d, and e (24, 38 
68, and 14 AU/10, 10, 7, and 10 Mr: lMarois et~aLl[2 008a. 
2010b), and Fomalhaut b (120 AU/< 2 M f lKalas et all 
20081) . While the sample of exoplanets is incomplete for 
a > 5 AU, indirect searches continue to hint that the 
semi-major axis distribution is at least flat, and possi- 
bly r ising, in dN/d\og{a) beyond 5 AU jC umming et al.1 
2008). The sample of six microle nsed planets be- 
yond the ice line supports this trend (jSumi et al.l 120101 : 
iGould et al.l [2010). Thus, a direct imaging search of 
outer solar system regions (5-50 AU), such as proposed 
here, would increase the total number of planets found 
relative to those in inner solar system orbits (< 5 AU). 
The goal of direct imaging is to assemble the first statis- 
tically significant sample of exoplanets that probes be- 
yond the reach of indirect searches and quantifies the 
abundance of solar systems like our own. 

2. OVERVIEW OF THE GPI INSTRUMENT 



The Gemini Planet Imager is configured to allow high 
contrast (> 10 7 ) imaging on angular scales of the diffrac- 
tion limit (~ 5 — 20X/D). High-contrast imaging with 
current AO systems is almost completely limited by 
quasi-static artifacts caused by slowly evolving wavefront 
errors. These originate from many sources, including in- 
adequately calibrated non-common-path errors that arise 
from the difference between the science and wavefront 
sensing paths, uncorrectable high spatial frequency er- 
rors on instrumental and t elescope optics, aliased high - 
frequency wavefront errors ([Povneer fc Macintoshl [2004) , 
chromatic errors, reflectivity variations, and Fresnel 
propagation effects. Although these can be partially 
attenuate d through techniqu es like angular differential 
imaging (Mar ois et al.l 1200a ). the errors evolve rapidly 
enough that ADI-like techniques do not operate well at 
small angles where it takes significant time for enough 
field rotation to accumulate. Unfortunately, these small 
angles are also the scales that correspond to our own 
solar system for most nearby young stars. GPI was de- 
signed from the beginning to minimize these quasi-static 
error sources and hence allow high-contrast imaging ap- 
proaching the photon noise limit. For more detail, see 
Macintosh etaLl (pool [2008h . 

The GPI system consists of five key subsystems: an 
AO system, a coronagraph, an interferometer, an integral 
field spectrograph, and a software system. 

The AO system makes fast measurements of the in- 
stantaneous wavefront and provides wavefront control 
via a tip-tilt stage and two deformable mirrors, one con- 
ventional piezoelectric and one high-order silicon micro- 
electro-mechanical-system (MEMS) device. The high- 
order AO system has 43 actuators across the diameter 
of the 7.8-m Gemini South primary mirror and operates 
at at an update rate of 1.5 kHz. The AO system uses 
a spatially-filtered Shack-Hartmann wavefront sensor to 
minimize aliasing an d have uniform response in varying 
atmospheric seeing (|Povneer fc Macintosh! I2004D . The 
anti-aliasing produces a characteristic square "dark hole" 
region, ~ 43A/Z? on a side. 

The coronagraph is an apodizer-pupil Lyot 
coronagraph (APLC) , and it controls diffraction 
(|Soummer et al.l [2006). This combines a carefully- 
designed input apodization with a focal-plane occulting 
stop (in GPI, a mirror with a central hole) and a 
Lyot pupil stop matched to the input telescope pupil 
size. The GPI APLC implementation is optimized for 
achromatic performance to improve multi-wavelength 
speckle suppression. Each waveband has a hard-edged 
occulter hole with a radius of 2.8 X/D, but that hole is 
also surrounded by a bright Airy ring, so high contrast 
is only fully practical at (~ 4-5 X/D). 

The interferometer provides low-temporal bandwidth 
precise and accurate measurements at the science wave- 
length of the time-averaged wavefront delivered to the 
coronagraph occulting spot. This information is used 
to remove slowly-evolving quasi-static errors due to flex- 
ure and changes in the response of the main visible-light 
wavefront sensor. 

The science instrument is an integral field spectrograph 
(IFS) that images in simultaneous multiple wavelength 
channels. The lenslet-based IFS has sampling of 0.014" 
per lenslet, with a square field of view of 2.8" on a side. 
Each spatial pixel is dispersed into a spectrum with reso- 
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lution R ~ 45. A single observation covers one of the Y, 
J, or H bands, or half of the K band, which is split in 
to K\ and K^. Raw detector images are reassembled by 
a data pipeline into 200 x 200 x 16 element data cubes. 
The IFS also includes a differential polarimetry mode, in- 
tended for characterization of circumstellar dust and not 
discussed in this paper. Finally, the software system co- 
ordinates communication between the other subsystems 
and the observatory software. 

GPI operates at a Cassegrain focus on the alt/az 
mount of the Gemini telescope in a fixed orientation with 
respect to the telescope, to increase instrumental stabil- 
ity and enable ADI post-processing. A typical science 
observation will consist of an hour long sequence of 30- 
60 second exposures. This is long enough that read noise 
and dark current are unimportant contributors to the 
noise budget, and short enough that image motion dur- 
ing the exposure is minimal. 

3. INSTRUMENTAL NOISE CALCULATIONS 
3.1. Simulations 

High-contrast imaging noise can be broken down into 
two main categories: speckle noise and photon noise. In- 
stantaneous monochromatic high-contrast images consist 
of a pattern of bright "speckles" surrounding the central 
core. These speckles have a size of ~ A/D, compara- 
ble to the image of a planet, and their random fluctua- 
tions are usually the main limitation in planet detection 
in existing high-contrast imaging instruments. In many 
cases (e.g. atmospheric turbulence wavefront errors), the 
speckle pattern will evolve rapidly, and after some char- 
acteristic speckle lifetime has elapsed the speckle noise 
will scale with exposure time as t 5 and in a long expo- 
sure produce a more uniform halo of light. In addition, 
any light present in the focal plane due to wavefront er- 
rors, whether speckled or uniform, will contribute Pois- 
son photon noise. 

In the GPI architecture, we consider the follow- 
ing sources of wavefront error and their correspond- 
ing speckle and/or photo n noise (|Povneer fc Macintosh! 
[20061 IMarois et al.ll2008bh : 

1. Residual atmospheric errors. Even operating at 
1.5 kHz, the motion of the turbulent atmosphere 
will cause changes in wavefront between measure- 
ment and the application of a correction. This is 
the dominant source of mid-to-low frequency wave- 
front errors. Here, low-frequency errors (< 3 cy- 
cles/pupil) refers to errors at spatial frequencies 
mostly blocked by the coronagraph occulter, and 
mid- frequency errors (3-22 cycles/pupil) are those 
transmitted by the occulter but that are within the 
controllable range of the deformable mirror. Errors 
at these frequencies contribute scattered light and 
photon noise near the star. Higher frequency spa- 
tial wavefront errors (> 22 cycles/pupil) are un- 
sensable and uncorrectable by the GPI AO system 
and contribute light outside the dark hole. A rep- 
resentative speckle produced by these atmospheric 
wavefront errors has a lifetime of a few tenths of a 
second and an amplitude of 60 nm. 

2. Wavefront sensor measurement noise. On dim- 
mer stars, the individual wavefront sensor measure- 



ments will include a significant noise component 
that will translate to random errors of position of 
the deformable mirror. In a Shack-Hartmann wave- 
front sensor, these errors are not spatially white, 
but have increasing power at low spatial frequen- 
cies. On dim stars, this is again a significant 
source of scattered light over the whole dark hole 
region. Speckles associated with this error source 
randomize at the closed-loop bandwidth of the sys- 
tem, with lifetimes of a few milliseconds, rapidly 
smoothing out. A typical value is 25 nm for a star 
with 7 — 6 mag. 

3. Residual non-common-path wavefront errors. Af- 
ter daytime calibration and closed-loop correction 
by the precision interferometric wavefront sensor, 
these are assumed to be ~ 5 nm of low-frequency 
error and ~ 1 nm of mid-frequency errors. The 
amount of light scattered by these small wave- 
front errors produces a negligible amount of photon 
noise, but the slowly-varying speckles they produce 
are a significant source of speckle noise. 

4. Reflectivity variations and amplitude errors. The 
GPI deformable mirrors are operated in phase 
conjugation mode, correcting the phase errors of 
the wavefront and producing a symmetric final 
PSF. As a result, any light scattered by changes 
in amplitude — for example, from reflectivity vari- 
ations on the Gemini primary mirror — is uncor- 
rectable. We assign a reflectivity variation to each 
individual optic in our system (typically 0.1%, ris- 
ing to 1% for the primary mirror and the MEMS 
DM, with a k~ 2 5 spatial frequency power spec- 
trum). The light scattered by these reflectivity er- 
rors again produces little photon noise but signifi- 
cant speckle noise. 

5. Fresnel propagation errors. Surface errors on an 
optic at an arbitrary location in the GPI opti- 
cal train will initially produce a pure phase error 
in the wavefront, but as the wavefront propagates 
these errors will mix bet ween phase and amplitude 
(|Shaklan fe Green|[2006l) . The finite size of GPI op- 
tics will also produce amplitude fluctuations near 
the edge of the beam. As noted above, GPI will 
only correct the component of these errors that is 
realized as phase at the deformable mirrors, not 
the amplitude component. To mitigate this, GPI 
optics are located away from focal planes and are 
of very high quality (~ 2 nm RMS wavefront error 
typically), but these are still a significant source of 
persistent and chromatic speckles. 

6. Atmospheric scintillation. Classical scintillation 
causes uncorrectable amplitude fluctuations in the 
telescope pupil similar to the Talbot propagation 
errors. Unlike the internal Talbot errors, those 
from the atmosphere are time varying. Simulations 
and calculations show that for our typical atmo- 
sphere profile this is a negligible contribution to the 
photon and speckle noise and hence it is ignored. 

In addition, normal astronomical noise sources such as 
sky background, detector readout, and dark current will 
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be present, though for typical GPI targets with H < 8 
mag. these are negligible. Telescope vibrations (wind- 
shake or mechanical vibrations of either GPI or the tele- 
scope secondary) are not included in this model, in part 
because the exact vibration environment of the Gem- 
ini telescope is unknown. The GPI tilt error budget 
allocates < 5 mas RMS to these effe cts; GPI employs 
an ad vanced tip/tilt control algorithm (jPovneer fc Veranl 
2010) that can cancel out resonant vibrations of the tele- 
scope, and the mechanical structure includes tuned-mass 
damping to cancel internal vibrations. These 5 mas have 
no significant effect on the coronagraph or on planet de- 
tectability. 

Modeling all these effects is challenging. Even with 
simple Fraunhoffer propagation, simulating the AO 
wavefront correction and propagation to the science fo- 
cal plane requires several CPU-hours per second of ex- 
posure time. Since GPI has been designed to mini- 
mize the quasi-static wavefront errors 3, 4, and 5 in 
the list above, their effects can only be seen after many 
minutes of integration. As a result, we treat the dy- 
namic and static wavefront errors as independent and 
treat them in two different simulations. The hrst simu- 
lates the dynamic behavior of the atmosphere and each 
component of the adaptive optics system to produce 
short-exposure point spread functions, using Fraunhof- 
fer propagation through the GPI optics and corona- 
graph. These were run for a standard Gemini C erro 
Pachon atmosphere (jTokovinin fc Travouiilonll2006l ) and 
for stars with brightness ra nging from / = 5-9 mag. 
dPovneer fc Macintosh! 120061) . The second, described in 
iMarois et al.l (|2008bf ). includes all the quasi-static wave- 
front error source and full Fresnel propagation through 
the GPI optical train, but has no atmospheric wavefront 
errors. Both simulations are run at each IFS wavelength 
channel within the H band. Photon and speckle noise 
are evaluated separately in each simulation. 

3.2. Post-processing 

Speckle noise can be attenuated through post- 
processing techniques such as ADI. However, this de- 
pends critically on the stability of the aberrations pro- 
ducing those speckles. We have assumed no ADI atten- 
uation whatsoever, making the pessimistic assumption 
that wavefront errors will evolve just fast enough to be 
unsubtractable over the timescales needed for ADI field 
rotation. Field rotation will still produce some averaging 
of the residual speckle noise as the planetary companions 
move through the speckle noise pattern, so we reduce the 
speckle noise by an amount equal to the square root of the 
number of X/D that field rotation will move the planet 
through in a 1-hour exposure. 

Since GPFs science instrument produces spectral data 
cubes, speckles can also be attenuated by scaling and 
subtracting different wavelength channels. We assume 
post-processing using a simple "double difference" algo- 
rithm using three different wavelengths. For the atmo- 
spheric speckles, which are almost purely phase errors 
and hence have well-behaved chromaticity, this attenu- 
ates speckle noise by more than 4 to 6 magnitudes at all 
radii, leaving the atmospheric speckle noise completely 
negligible. For the quasi-static speckles, especially the 
Fresnel effects, the speckles are more chromatic and they 
are attenuated by only 2-3 magnitudes within the dark 



hole region. Since these speckles are weak to begin with, 
the residual noise is lower than the photon noise for all 
but the brightest stars. 

This double-difference subtraction will also attenuate 
the planetary signal, unless the planetary signal con- 
tains deep molecular absorption features (e.g. methane), 
or the wavelength scaling of the speckles cause them 
to move by more than 2X/D through a single spectral 
band, i.e. at radii greater than 10X/D for a 20% band- 
pass. Although not all planets will show methane ab- 
sorption, the faintest (and hence most difficult to de- 
tect) planets are expected to. The HR 8799 planetary 
companions do not show methane absorption but are so 
bright — contrast ~ 10 5 relative to their F-star host — 
that they would be easily detected by GPI even with- 
out spectral differencing. Hence the use of the spectral- 
difference contrast curves is appropriate for the Monte 
Carlo modeling described here. In actual operation, GPI 
will likely use a more sophisticated approach to PSF sub- 
traction with wa velength and time, like ly based on the 
LOCI algori thm dLafreniere et~aLl l2007aft and the SOSIE 
framework (jMarois et al.l I2010af) . However, this would 
not significantly change the results here, since in these 
simulations even with simple spectral differencing GPI is 
photon-noise rather than speckle-noise limited for stars 
brighter than 1 = 6 mag. 

3.3. Final instrumental contrast predictions 

The results of the instrumental noise calculations are 
tabulated as a function of stellar brightness and angular 
separation for each noise source. These tabulated values 
are used to generate the total number of photons con- 
tributed by all sources of noise for some exposure time. 
A planet may then be considered detectable if the ex- 
pected number of photons detected from the planet is 
five times the number of photons contributed by noise. 
The contrast curves in Figure [T] summarize this result, 
showing the ratio of five times the overall noise to the 
number of photons detected from the host star as a func- 
tion of angular separation. The instrument is expected 
to reach contrasts better than 10~ 7 over an hour long ex- 
posure for targets with / ~ 5 mag. For brighter targets, 
the performance remains the same, while the expected 
contrast declines for dimmer targets. For stars beyond 
/ ~ 10-11 mag., the system performance becomes com- 
parable to existing instruments. For that reason, stars 
brighter than / ~ 10 mag. will make up the majority of 
GPI targets. 

3.4. Matched filters 

Models of the spectral properties of planets (these 
models will be discussed in more detail in the next sec- 
tion) predict young planets will show many spectral fea- 
tures, rather then being uniformly bright across a wave- 
length band. For this reason, the integral field spectro- 
graph, which provides 16 spectral channels across each 
wavelength band, will be important to the success of GPI 
at detecting planets. A planet may be bright enough in 
some spectral channels to be above the contrast limit, 
while its flux when summ ed over the en t ire ba nd is not. 
The model spectrum from lFortnev et ahl ([2008D in Figure 
[2] is an example of such a case. The flux from this model 
planet is concentrated in a narrow range of wavelengths, 
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Fig. 1. — Contrast curves (5<r) in the H band as a function of 
angular separation for five different star I magnitudes, assuming 
one hour on the target and an AO star. Photon shot noise domi- 
nates for stars with / > 6 mag., while static speckles dominate for 
stars with I < 6 mag. Visit http://planetimager.org for a tabular 
version of these data. 




Wavelength (fJ.m) 

Fig. 2. — A smoothed model spectrum for an atmosphere with 
T c ff = 400 K and log(g) = 3.66 is shown as a solid black line, with 
black points indicating where the spectrum is sampled, and system 
noise in each channel as solid black bars. This is representative of a 
4 Mj planet around a 500 Myr Sun-like star 40 parsecs from Earth. 
The dashed black line shows terrestrial atmospheric transmission. 



while the rest of the spectral channels contribute primar- 
ily noise. When looking at the spectrum, even without 
the model spectrum to guide the eye, it looks like a clear 
detection. Yet the overall signal to noise ratio across 
the entire H band is less than 5, meaning a sum over all 
spectral channels would yield a non-detection. 

Analysis of real data will involve much more than eval- 
uating a planet's signal to noise ratio and declaring a de- 
tection or non-detection. For the purposes of simulating 
tens of thousands of planets around thousands of stars 
though, this is a practical approach. Yet the spectral fea- 
tures in planetary atmospheres suggest that generating a 
planet's signal to noise ratio by taking a simple sum of all 
spectral channels will underestimate GPI's ability to de- 
tect that planet. To explore how much the detection rate 
may be affected, we used model atmospheres to generate 
matched filters. This was done by degrading the resolu- 
tion of a model atmosphere at a given surface gravity and 
temperature to match that of GPI (A/AA ~ 45 in the H 
band), and then weighting each of the spectral channels 
according to the model predicted flux. These weighted 
channels are then summed to find the total flux from the 
planet, and the total noise. 

Using matched filters significantly improves survey per- 
formance. Relative to a filter that gives equal weight to 
all spectral channels, a matched filter corresponding to 
an atmospher e at 400 K and surfac e gravity log (3) ^3.6 
(in cgs) from iFortnev et al.l (|2008t ) provides a factor of 
3 improvement of the detection rate. This assumes that 
the pl anet atmospheres are as described bv lFortnev et al.l 
(2008) as well, so this represent s the ideal situation 
Using the model atmospheres of iBurrows et al.1 ()2003[ ) 
to generate a matched filter instead, the result changes 
marginally. However, given uncertainties in planet at- 
mospheres, particularly cooler gas giants, all subsequent 
simulations results reported here do not use matched fil- 
ters. As a consequence, estimated planet detection rates 
may be somewhat pessimistic. 

4. PLANET MODELS 

We consider planets without internal energy sources 
(i.e., deuterium fusion) that are young and sufficiently 



far from their primary star that we can ignore the stel- 
lar radiation field when considering cooling and contrac- 
tion, i.e., T off >> (i?»/a) 1/2 r*, or T cff 4 >> (180iT) 4 for 
a planet orbiting a solar type star at a — 5 AU. 

The luminosity of a planet is given by its radius and 
effective temperature, but the detectability at specific 
wavelengths also depends on the nature of the dominant 
opacity sources in the photosphere. The photospheric 
chemical composition at a given effective temperature 
is sensitive to surface gravity, g — GM / R 2 . Both T e g 
and log(<?) are depend upon the age and the mass of the 
planet. To understand the detectability of a planet then 
involves converting a planet age and mass to a luminosity, 
effective temperature, radius, and surface gravity. 

To establish this mapping, we compare two evolution- 
ary models. One is the so called "hot start" model of 
IBurrows et al.l (119971) . and the other is the the core ac- 
cretion model of iMarley et al.l ((2007). The most signifi- 
cant difference between the two is summarized in Figure 
02 which shows effective temperature and surface gravity 
for different mass and age contours. While the expected 
effective temperatures and surface gravities of the two 
models come in to rough agreement for older planets, 
planets that form via core accretion are expected to be 
significantly cooler post-formation than planets formed 
via a hot start. 

4.1. Hot start models 

The label "hot start" was used bv lMarlev et al.l (|2007l ) 
to describe models of the atmospheres and evolution of 
extrasolar giant planets that ignored the different forma- 
tion mechanisms of brown dwarfs and gas giants. These 
models assumed that an object's mass, age, and com- 
position dominated its spectral characteristics, and as- 
sumed that gas giant planets, like brown dwarfs, form 
like stars. Objects formed in this manner possess a fully 
adiabatic interior with high specific entropy, which cor- 
responds to a high internal temperature. These models 
focus on treating the atmospheres of gas giant planets 
in order to predict the emergent flux as a function of 
temperature and surface gravity. 
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Fig. 3. — Surface gravity plotted against effective temperature for 
selected mass and age contours for two planet for mation models 
where contours are labeled in Mj and log(t). ThelBurrows et~al 



fl2003l ) hot start model is shown in black, and the IMarTeT^t*^T 
H2007T ) core accretion model is in cyan. The H2O condensation 
line is shown in red, below which planet atmosphere calculations 
become more uncertain. The position of Jupiter, at T e jf i=a 125 K, 
log(g) £S 3.4, and log(i) ss 9.66, is shown as a magenta dot. 



The hot start model used here is that of lBurrows et al.l 
(fl997h . which spans an effective temperature range from 
approximately 100 K to 1000 K, covering the realm be- 
tween the known Jovian planets and the known T dwarfs. 
The models can be conveniently approximated as power 
laws, found by performing a regression over a restricted 
range of masses and ages of interest: 1 < M/Mj < 12 
and 0.01 < t/Gyr < 2. The resultant expressions are 
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These expressions are comparable t o previous se miempir- 
ical estimates of these dependences (Black 1980), and the 
luminosity expression is close to that of cooling at con- 
stant heat capacity (L ~ Af 2 i~ 4 / 3 ). The resultant radii 
and effective temperatures have rms errors of 2% and 
4% respectively. The corresponding error in bolometric 
magnitudes is 0.15 mag. We use th e radiative-co n vectiy e 
equilibrium atmosp here model o f Marlc vet al.l (|1996f) . 
further described in lBurrows et al.l (11997); the most re - 
cent tabulations are provided in iBurrows et al.1 ([2003). 
The model has been updated to self - consis tently include 
both alkali op acities IBurrows et al. (I2000D and precipi- 
tating clouds (|Ackerman fc Marlevll200l( )7 

A major source of uncertainty in planet atmosphere 
models is the treatment of clouds. As with the known T 
and L dwarfs, absorption by water vapor dominates the 
spectra of the cooler brown dwarfs. These features gen- 
erally deepen with increasing age and decreasing mass. 
This trend is in part due to the increase with decreas- 
ing gravity of the column depth of water above the pho- 
tosphere. At effective temperatures below 400-500 K, 
water vapor condenses in planetary atmospheres. The 



appearance of water-ice clouds constitutes a major un- 
certainly separating the known T dwarfs from the giant 
planets, and is denoted in Figure [3l When water con- 
denses, water vapor is depleted above the cloud tops 
causing a decrease at altitude in the gas-phase abun- 
dance of water. Water clouds form in the atmosphere 
of an isolated 1 Mj object within 100 Myr, and within 
2 Gyr they form in the atmosphere of a 12 Mj object. 
However, at supersaturations of 1% and for particle sizes 
above 10 fim, such clouds (and the corresponding wa- 
ter vapor depletions above them) only marginally affect 
the calculated emergent spectra. For wavelengths long- 
ward of 1 /zm, the cloudy spectra differ from the no-cloud 
spectra by at most a few tens of percent. Below effec- 
tive temperatures of 160 K, NH3 clouds form. This is 
likely well below the effective temperature we can hope 
to detect. 

We expect that planets detected in reflected light will 
comprise only a negligible fraction of our sample, and 
therefore we have not attempted a thorough treatment 
of this problem. Rather, we assume a Bond albedo of 
0.4. This assumption would gives rise to errors if the ob- 
serving band were coincident with a strong CH4 band or 
when the target it is a so-called Class III or "clear" extra- 
solar giant planet, so named because they are expected 
to be too hot (T e g> 350 K) to c ontain any principal con- 
densates (jSudarskv et al.ll2000ft . 

4.2. Core accretion 

The standard theory for the formation of gas giant 
plane ts is the core accretion model (e.g., iPollack et al.l 
1996), which begins with dust particles colliding and ag- 
glomerating within a protoplanetary disk to form icy 
and rocky planetary cores. If the core becomes mas- 
sive enough while gas remains in the disk, it can grow 
by gravitational accretion of this gas. Gas giants ac- 
crete most of the gas within their tidal reach, filling the 
Hill sphere around them with a hot, extended, gaseous 
envelope. Further accretion is slowed by the dwindling 
supply of local raw materials and by the extended en- 
velope, leading to growth times of 5-10 Myr. The pre- 
dicted planet formation time is uncomfortably long com- 
pared to the observed ~ 3 Myr lifetime of protoplan- 
etary disks. Two factors may alleviate this time-scale 
problem: 1) inward migration can bring giant planets 
to fresh, gas-rich regions of the disk; 2) dust sedimen- 
tation may reduce atmospheric opacity, which leads to 
more rapid escape o f accretion lum i nosity and shrink- 
age of the envelope. iHubickvi et al.l (|2005| ) have shown 
that redu cing the grain opa city to a level observed in 
L dwarfs (iMarlev et al.|[2002l ) reduces the planet growth 
time scale to ~ 1 Myr, well within the lifetime of proto- 
planetary disks. Once accretion stops, the planet enters 
the isolation stage and the planet contracts and cools at 
constant mass. 

The configuration of the protoplanet at the end of run- 
away gas accretion represents the i nitial conditions for 
subs equent cooling, a nd co ntraction. IMarlev et al.l (|2007l ) 
and iFortnev et al.l (|2008l ) have conducted preliminary 
calculations that describe the cooling and contraction of 
a young planet as it emerges from its parent disc. The im- 
plication of the IMarlev et al.l (|2007l) results is that giant 
planets formed by core accretion are less luminous post- 
accretion than had been previously expected because sig- 



nificant energy is radiated during the formation process. 
The fully formed planet has a sm aller radius at young 
ages than hot start models predict (|Burrows et al.ll 19971 : 
iBaraffe et al.1 120031 ) . leading to a lower post-formation 
luminosity. There are two significant observational con- 
sequences: 1) there is a period of very high luminosity 
(~ 10~ 2 L ) which lasts ~ 40,000 yr; 2) the initial con- 
ditions for subsequent evolution are not "forgotten" for a 
time of order the Kelvin-Helmholtz timescale, which lasts 
tens of millions of years. These factors imply that obser- 
vations of class II (0.5-3 Myr) and III (~ 10-100 Myr) 
young stellar objects afford the opportunity to probe the 
planet formation events. The run-away accretion spike 
is likely to be broader and fainter than in these idealized 
calculations because of gradual accretion across the gap 
that the protoplanet forms, and the probability of wit- 
nessing this event in a typical T Tauri star may be much 
larger than a few percent 

The approximations used to compute the luminosity 
history of a planet formed by core accretion are similar to 
those adopted in early studies of protostellar formation — 
the runaway gas accretion phase is described using the 
formalism of iStahler et al.l (|1980f) , where matter falling 
onto the central object passes through a 1-d, optically 
thick shock. The 1-d geometry forces all the matter to 
pass through this shock, whereas in nature the accretion 
is at least 2-d and more realistically 3-d. Comparison 
with star formation suggests that some accretion may 
occur on viscous time scales in a disk r ather than on 
the fa st dynamical timescale depicted in iFortnev et al.l 
(2008). Ultimately, the virial theorem must be satisfied, 
and half the grav i tation al potential energy is radiated: 



the IMarlev et al.l (12007) calculations and the ho t start 



models (Burrows et al 



19971: [Baraffc c t~aT1[200l repre- 



sent limiting cases in the contraction histo ry. 

Des pite the preliminary nature of the IMarlev et al.l 
(2007) results, one message is clear: the luminosity 
of young exoplanets encodes key information about 
how they were formed. Simulations suggest that the 
timescales for relaxation are longer for more massive 
planets, and are at least approximately related to the 
Kelvin-Helmholtz timescale, with 



f2 



GM 
RL ' 



(4) 



where G is the gravitational constant, M is the mass of 
the planet, R is the radius of the planet, and L is the lu- 
minosity of the planet. This timescale ranges from tens 
to hundreds of millions of years for gas giants, meaning 
observations of stars with ages in this range will yield 
information about planet formation. Young targets are 
also the most promising targets because a brighter planet 
is more easily detected, but planets formed via core accre- 
tion may be cooler at young ages than hot start planets. 
If planets are significantly cooler at young ages than hot 
start models assumed, the detection rate of young gas 
giants will be lower than expected. 

5. OBSERVATION PLANNING 

5.1. Planet populations & Monte Carlo simulations 

We use the results of the sensitivity calculations de- 
scribed in Section [3] and the planet formation and atmo- 
sphere models described in Section [4] to perform Monte 



Carlo simulations characterizing the resu lts of different 
surveys using the Gemini Planet Imager ()Graham et al.l 
l200j TGraham 20071 [Graham et al.l 120071) . The sensitiv- 
ity of the instrument to a planet depends upon the an- 
gular separation between the host star and planet and 
the brightness of the host star. The spectrum of a planet 
depends on its effective temperature and surface grav- 
ity, which themselves are functions of the mass and age 
of the planet. Specifying the mass, age, and orbital el- 
ements of a planet and the brightness and distance of 
its host star is then sufficient to determine whether that 
planet is detectable according to the model used. By gen- 
erating distributions of the orbital elements and mass of 
exoplanets, and assuming that planet and host star for- 
mation are coeval, we may then simulate the likelihood 
of detecting a planet around some star. 

Lacking an accepted paradigmatic theoretical model 
for the distributions of these parameters, we rely on the 
distributions at semi-major axes smaller than ~ 5 AU 
found via radial velocity surveys, and extrapolate these 
distributions to larger semi-major axes. The applica- 
bility of the mass and eccentricity distributions found 
for close-in companions at large semi-major axes is un- 
clear, as is the extension of the semi-major axis distribu- 
tion to larger semi-major axes. Nevertheless, as a base- 
line we apply recent results that indicate that the planet 
mass and semi-major axis each follow a power law dis- 
tribution, with dN/dM oc M -131 and dN/da oc a~°- 61 
(jCumming et al.ll2008| ). For the mass, we adopt a lower 
cutoff of 0.5 Mj, below which a planet is so dim as to 
make direct imaging unlikely, and below which radial ve- 
locity surveys are incomplete. We take as an upper cut- 
off 12 Mj, reflecting the transition from planet to brown 
dwarf at the deuterium fusion limit, thereby ignoring any 
details of how the planet formed. A histogram of planet 
masses sampled from this distribution is shown in the 
upper left panel of Figure 2] 

Similarly, we adopt cutoffs for the the semi-major axis 
distribution of 0.1 AU and 75 AU. Radial velocity sur- 
veys show a rising power law for the semi- major axis dis- 
tribution in log space out to ~ 5 AU, and iNielsen et al.1 
IpOlCh placed an upper limit of 75 AU on the s emi-m ajor 
axis distribution quoted in iCumming et al.l (|2008l ). A 
lower cutoff may be justified, but already detected exo- 
planct systems, such as Fomalhaut, HR 8799, and 1RXS 
J 160929 indicate that some planets form at large semi- 
major axes. Such a cutoff i s also consistent with the 
typical sizes of T Tauri disks (jlsella et al.ll2009h . A sam- 
pled population from this distribution is shown in the 
upper right panel of Figure [¥] 

The eccentricity distribution is also drawn from radial 
velocity observations. Over a decade ago, an ecc entricity 
distrib ution of dN/de oc e -0 5 was proposed bv iHeacoxl 
(fl999T ). but to our knowledge, no updated eccentricity 
distribution has been published. We found that a sim- 
pler linear fit to eccentricities adequately describes obser- 
vations to date, and adopted that, as results are insen- 
sitive to the eccentricity distribution anyway. The rest 
of the orbital elements are easily obtained. The mean 
anomaly and the argument of perihelion are both drawn 
from a uniform distribution from to 2ir. The incli- 
nation of the orbit relative to the plane of the sky is 
drawn from a distribution of randomly oriented orbits, 
i.e. dN/di = \ sin(i). 
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Fig. 4. — Distributions of planet properties in a typical survey. The mass and semi-major axis distributions are drawn from extrapolations 
of distributions consistent with radial velocity detected planets. The angular separation distribution assumes a population of host stars 
distributed randomly in a volume within 100 parsecs of the Sun. The effective temperature distribution assumes a distribution of stellar 
ages based on a sim ple galactic disk star formation rate, and that planet mass and age map to effective temperature as described in 
IMarlev et all J2007T) . 



With these distributions, we use a pseudo random 
number generator and the rejection method to generate 
a population of exoplanets, which may be placed around 
some star. The distribution of masses is then converted 
to distributions of surface gravity, radius, and tempera- 
ture according to the adopted formation model, and the 
temperature and surface gravity then give the emergent 
spectrum of the planet, as described in Section 2J From 
this, an expected flux at the telescope in each spectral 
channel is calculated. Per the discussion in Section 13.41 
these channels may be summed with our without weights 
to get a total flux across the band, but in all results here 
channels are summed without weights. Finally, the dis- 
tributions of orbital parameters yield a distribution in 
angular separation. For a given star, we then have rel- 
ative brightnesses as a function of angular separation, 
which can be compared with the expected performance of 
GPL The probability of detecting a planet around a star, 
assuming a single planet with mass 0.5Mj < M < 12Mj, 
is then directly obtained from the percentage of relative 
brightnesses across the entire wavelength band that are 
above the threshold of detectability, which is taken to be 
a signal to noise ratio of 5. 

Known exoplanets suggest that these distributions 
should depend on the properties of their host star. The 
observed metallicity dependence for short period plan- 



ets has already been discussed in Section [TJ but will 
not be included in any simulations due to uncertainty 
in its applicability at large semi-major axes. Doppler 
surveys show that the fraction of stars with planets in- 
crease s with stellar mass over the range from M to A 
stars (| Johnson et al.ll2008l ) . Microlensing results for the 
distribution of planets beyond the ice line are in terms 
of the mass ratio distribution dN/dq, where q — M p /M s , 
rather than a simple distribution of masses dN/dM. The 
success at finding planetary mass companions around A 
stars with other direct imaging efforts also supports the 
importance of host star mass, even a t large semi - major 
axes (e.g., Fomalhaut and HR 8799). iGorti et all (2009) 
predict that the disk lifetime is relatively constant for 
M < 2M@, which is also consistent with the notion that 
Jupiter mass planets are more common around A stars 
than they are around M stars. To explore the impact 
of host star mass on direct imaging surveys, we tested 
both a mass ratio distribution dN / dq and a simple mass 
distribution dN I dM. T he mass power law index given in 
iCumming et al.l ()2008f ) was also used for the mass ratio 
distribution, so that dN/dq oc g -1- . The upper and 
lower cutoffs for q are set such that the lower mass cut- 
off is 0.5 Mj and the highest mass is 12 Mj, for the 
same reasons those cutoffs are adopted for the dN/dM 
distribution. This is unsubstantiated quantitatively, but 
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seems reasonable qualitatively. For the rest of the paper, 
we note when results are for a mass ratio distribution. 

Due to incompleteness, Doppler surveys have not 
yielded a n absolute planet hos ting probability. Nev- 
ertheless, iCumming et al.1 (|2008| ) estimate that 17-20% 
of FGK stars host planets with M > 0.3M,/ and with 
a < 20 AU. Surveys using Spitzer IRAC and MIPS find 
debris dis ks around 15% of solar -type stars younger than 
300 Mvr lCarpenter et all (f2009h and around 33% of A- 
type stars ISu et al. I (|2006t) younger than 850 Gyr, sug- 
gesting a lower limit of this order for the probability of 
hosting a planet. We generally sidestep the uncertainty 
in planet hosting probabality by reporting the detection 
rate in terms of the ratio of the number of planets de- 
tected to the total number of planets. When reporting 
the total number of planets detected around a sample of 
a fixed number of stars, we adopt the simple assumption 
that all stars host a single planet more massive than 0.5 
Mj in the range 0.5-70 AU. 

5.2. Simulated surveys 

Simulating a survey is a straightforward extension of 
simulating the likelihood of detecting a planet around 
a single star. Two of the authors independently wrote 
code to do this, and used slightly different approaches. 
In one approach, a sample population of planet proper- 
ties was drawn from the distributions described in the 
previous section and then placed it around each star in 
some sample of stars. In the other approach, a new sam- 
ple population of planet properties was created for each 
star. The results of these two approaches were in ex- 
cellent agreement, providing confirmation that neither 
code contained major errors, and the choice of using the 
same planet property distribution for each star or dif- 
ferent ones was unimportant to the results for simulated 
surveys with a large number of planets. 

The stellar sample used in a simulation may be either a 
list of real stars, or a Monte Carlo population, where stel- 
lar masses are generated from a Kroupa IMF, stellar ages 
are based on a simple galactic disk turn-on and turn-off, 
and stellar distances are generated to matc h the observed 
stellar density in the so lar neighborhood (|Kroupall2001l : 
iMihalas fc Binnevlll981[ ). In either case, the same proce- 
dure for simulating the likelihood of detecting a planet 
around a single star may be used for a number of stars. 
For different stars, the primary factors that change are 
the angular separation on the sky of star and planet, the 
brightness of the star, and the brightness of a planet of 
given mass, as planet brightness depends upon the age 
of the system. A record is kept of the planets that are 
detectable around each star. 

Performing a simulation over a sample of stars provides 
a total detection rate and a distribution of properties for 
detected planets. These may be used to compare sur- 
veys emphasizing different target selection, which may 
then inform where target identification efforts should be 
focused. With this in mind, we simulated a population 
of stars within 100 parsecs, and then performed three 
surveys for samples where targets were selected based on 
different sets of criteria. These criteria were age, dis- 
tance, and spectral type (or mass), and the results are 
shown in Figure [5] To provide a fair comparison between 
the surveys, the criteria adopted for each survey were 
chosen to yield roughly the same number of available 



targets in each survey. The distributions plotted assume 
complete information about the planet detected, which 
will obviously not be the case for real planets; while the 
semi-major axis and effective temperature may be de- 
rived from observations directly, determining the mass of 
a planet from direct observations will generally require 
planet formation and atmosphere models, or dynamical 
constraints for multi-planet systems or planets in disks. 

The results in Figure [5] show that observing young 
stars that are even reasonably near the Sun may be ex- 
pected to yield the highest detection rates, as well as the 
most information about planets at the lower end of the 
mass distribution. The primary drawback in observing 
the youngest stars is that it requires observing stars at 
greater distances. The inner edge of the "dark hole" re- 
gion is roughly 0.2", so only planets outside of 20 AU may 
be observed for a star that is 100 parsecs distant. If gas 
giants are very rare at semi-major axes greater than 20 
AU or so, selecting only young stars would lead to a sig- 
nificantly diminished detection rate. Selecting for stars 
that are somewhat older and nearer is not as successful 
in detecting planets for the semi-major axis distribution 
that we have assumed, as the gas giants with masses 
~ Mj and ages ~ Gyr are too dim to detect in an hour 
observation in most cases. Such a survey choice does 
have the advantage that the majority of detected plan- 
ets are at semi-major axes < 20 AU, and the existence 
of gas giants at 5-20 AU is somewhat better established 
than the existence of gas giants from 30-70 AU. Present 
efforts to identify targets for GPI are focused primarily 
on young stars, but very nearby stars that also happen 
to be only moderately young are promising targets as 
well. The desirability of early type stars depends on 
how the planet mass distribution depends upon stellar 
mass. Under the assumption that the planet mass dis- 
tribution is independent of stellar mass, later type stars 
have the highest planet detection rate, since for the same 
age and planet mass, the ratio of planet brightness to star 
brightness will be higher for dimmer stars. Planet mass 
very likely does depend on host star mass, and making 
the simple assumption that planet mass scales linearly 
with host star mass, planet detection rates are highest 
for early type stars. This result suggests a continued 
emphasis on identifying early type stars will be valuable 
to a successful direct imaging survey. 

Another means of survey comparision is the com- 
pleteness diagram, which encapsulates information about 
both the rate of detection and the properties of the plan- 
ets detected, showing for a grid of planet masses and 
semi-major axes the fraction of planets that are detected. 
Two samples are shown in the completeness diagram in 
Figure El with one sample limited to stars younger than 
100 Myr, and the other limited to stars younger than 1 
Gyr and nearer than 50 pc. The young star sample was 
taken from catalogs assembled by I. Song and M. Bessell 
(private communication). Recognizing the importance of 
a large target set of young, nearby stars (defined as age 
< 100 Myr and distance < 75 pc), Song and Bessell initi- 
ated a large-scale program in 2007 to identify and charac- 
terize young stars. There are approximately 200 known 
nearby yo ung stars (< 75 pc; < 100 Myr) in the litera - 
turc (e.g- lZuckerman fc Son g 2004; Torres et al.ll200l . 
To find new candidate young stars, Song and Bessell se- 
lected about 3000 bright (/ < 9 mag.) Tycho-2 stars from 
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Fig. 5. — Each row of panels plots the intrinsic mass, semi-major axis, and temperature distributions of detected planets from three 
different surveys. The target selection criteria for the three surveys in each row are given by the matching color text in the left column, 
and the total detection rate for the survey is given in the middle column. The black line is the same in all rows, showing a volume limited 
selection. Within each row, the cyan and magenta lines represent the same target selection criteria, but are distinguished by whether 
planet mass depends upon host mass. For cyan lines, planet mass is independent of host mass, while planet mass scales proprotionally to 
host mass for the magenta lines. Cyan and magenta lines are for a survey emphasizing massive stars in the top row, a survey emphasizing 
nearby stars in the middle row, and a survey emphasizing young stars in the bottom row. 



ROSAT catalogs with enhanced X-ray emission and po- 
tential young star kinematics and obtained optical echelle 
spectra of about 2000 high priority targets with the 2.3- 
m at Siding Spring. These spectra were used to extract 
age indicators including Li 6708 A, Ca II HK, Ha, and 
v sm(i ) and estimated ages following I Zuckerman &: Songl 
(2004), yielding about 200 additional young stars. When 
combined with an unpublished list of nearby young stars 
(400; Song, Zuckerman, and Bessell since 2000) and cur- 
rently known members of nearby young stellar groups 
(400), there are approximately 1000 distinct solar- type 
stars that are 100 Myr old or younger within 75 pc. There 
are also more than 1000 adolescent stars (100-500Myr.) 
Age uncertainties in the age-dating methods used are age 
dependent (smaller for younger ages) and ±5 Myr for the 
youngest stars (~10 Myr old) and ±300 Myr for the old- 
est stars (500 Myr). The sample used here is limited to 
roughly 600 stars accessible to the Gemini South tele- 
scope, and has a median age of 50 Myr and median dis- 
tance of 50 pc. The mean age and distance of the sample 
are both 20% larger than the median. The other sam- 
ple consists of stars from the Geneva-Copenhagen survey 



(jHolmberg et al.ll2009T ) . which was chosen because it is a 
publically available catalog of stars with estimated ages 
and distances. From the Geneva-Copenhagen catalog, 
we selected stars accessible to Gemini South that are 
younger than 1 Gyr and nearer to the Sun than 50 pc, of 
which there are roughly 200. The diagram shows that a 
survey of moderately young and nearby stars will probe 
a reasonable fraction of mass and semi-major axis space, 
but only by surveying a very young sample of stars is 
there a significant probability of imaging gas giants with 
M ~ Mj. There is, however, a reasonable likelihood of 
detecting the more massive gas giants around adolescent 
stars, and detection of the most massive planets will be 
limited primarily by the fraction the orbit during which 
the planet is in the "dark hole" observing region of GPL 

5.3. Target Ordering 

In simulated surveys, a planet detection probability is 
recorded for each star. These detection probabilities are 
some function of the star's age, distance, and spectral 
type. We find this function can be approximated reason- 
ably well with a product of power laws when considering 
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Fig. 6. — Completeness diagram comparing two surveys that 
assum e planets evolve according to the models of IMarlev et al.l 
l|2007fl . The cyan represents detection rates for a survey of stars 
within 50 parsecs and younger than 1 Gyr target stars from 
the Geneva-Copenhagen Survey (Holmbcrg ct al. 2009), while the 
black lines are for a survey of young stars (age < 100 Myr) from 
a target list being compiled for GPI (Song and Bcssell, private 
communication). The colors of the text correspond to the same 
surveys, and show the number of stars surveyed N* and the num- 
ber of detected planets N p , assuming one planet per star. 
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Fig. 7. — The expected number of detected planets is shown for 
two observing strategies. The dashed black line indicates targets 
chosen randomly from a sample of stars younger than 1 Gyr and 
within 70 parsecs. The solid black line shows the results when the 
targets are ordered according to Equation [5] 

a ~ —5,(3 ~ —1.5,7 ~ 1- Likewise, changing the slope 
on the semi-major axis distribution unsurprisingly im- 
pacts the importance of stellar distance on a survey. 



only stars younger than 2 Gyr. For an older sample, the 
large number of stars with effectively zero probability of 
observing a planet can skew the resulting fit. Denoting 
the star's age by t, distance d, and mass M, the detection 
probability p is 

with the parameters scaled to appropriate values to ease 
fitting, and the logarithm of the age being used as the 
ages for young stars are known with precision at the or- 
der of magnitude level. Performing this regression over 
multiple trials, and assu ming that planet plan ets are de- 
scribed by the models of IMarlev et al.l (|2007l ) and that 
mass is independent of stellar mass, we find a ~ —5, 
(3 ~ —1.5, 7 ~ —1.5. While the best fit values are sensi- 
tive to the average age of stars in a sample, they empha- 
size priority in the same general way: the log of the age 
is the most important parameter, and the distance and 
mass of the star are of roughly equal importance. 

Once generated, this detection probability function can 
be used to order target selection before GPI begins ob- 
servations, as well as for ordering targets in new simu- 
lated surveys. This application is demonstrated in Figure 
[3 which shows the results of two surveys of 1000 stars 
limited in age to be younger than 1 Gyr. If all stars 
in the sample are observed, the end result is the same, 
but when ordering the targets, approximately 2/3 of de- 
tectable planets are found in the first 1/3 of the sample. 

Fitting for the target ordering parameters also gives 
an at-a-glance sense of what is important in a survey. 
As such, it is an interesting way to test the impact 
of assumptions about the planet population. For in- 
stance, when the distribution dN / dq is used rather than 
dN/dM for generating planet masses, detection proba- 
bility increases with host star mass rather than decreas- 
ing. The other parameters stay roughly the same, with 



6. EXPECTED GPI PERFORMANCE 
6.1. Detection rate and distributions 

With 150 nights of telescope time, GPI could perform 
initial 1 hour observations and follow-up for roughly 1000 
stars. The results of such a survey depend strongly on 
target selection. The predictions of the formation model 
also strongly influence an exoplanet survey, though this 
effect decreases as the mean age of targets increases and 
model predictions converge. An age limited survey of 
stars younger than 1 Gyr w ithin 80 parsecs wil l have a 
detection rate of 4% for the IMarlev et all (120071) model, 
where as planets that evolve according to iBurrows et al I 
(|2003f) will be detected around 12% of stars. For a sur- 
vey of the stars younger than 100 Myr being compiled by 
Song (2010, private communication) for GPI, and with 
the important assumption that the semi-major axis dis- 
tribution observed for radial velocity exoplanets contin- 
ues out to ~ 70 AU, these dete c tion r ates i ncrease to 13% 
and 2 1% for the IMarlev et al.l (|2007l ) and IBurrows et al.l 
(2003) models, respectively. For a volume limited survey 
of very nearby stars, the detection rate drops to 1-2% 
for both models. 

While the existence of a significant population of plan- 
ets at large semi-major axes (30-70 AU) will certainly 
improve the detection rate of a GPI survey, GPI may still 
be successful so long as gas giants exist at more moder- 
ate semi-major axes. GPI is capable of detecting planets 
separated by less than 0.2" from their hosts. This is il- 
lustrated by Figure [SJ which shows the distribution of 
angular separation for planets detected by GPI, as well 
as the intrinsic angular separation distribution, given the 
assumed semi-major distribution and distribution of stel- 
lar distances. The majority of detected planets will come 
from the inner half of the GPI "dark hole." For the in- 
nermost angular separation bin, the average semi-major 
axis of detected planets is 18 AU. These small separa- 
tions are where GPI will be most valuable, probing re- 
gions of semi-major axis space inaccessible to previous 
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Fig. 8. — The distribution of angular separation for GPI detected 
planets is shown by the solid black line, while the distribution for 
generated planets is shown by the dashed black line. The average 
semi-major axis of planets within each bin is printed at the top of 
the bin. 



imaging surveys. Figure also indicates that regardless 
of the number of planets at large semi-major axes, GPI 
will probe a region of semi-major axis space which has 
been unexplored by previous radial velocity surveys and 
direct imaging surveys. 

Likewise, GPI is capable of detecting a wide range of 
planet masses. As shown in Figure [5j observations of 
young stars may be expected to yield detections of plan- 
ets with mass M ~ Mj in reasonable numbers. For plan- 
ets with mass M ~ 8Mj and above and ages less than 
1 Gyr, the detection rate is limited primarily by frac- 
tion of an orbit that planets spend in the "dark hole" 
region. This remains true for objects beyond the deu- 
terium burning demarcation between planets and brown 
dwarfs. Imagining the planetary mass distribution ex- 
tending to the realm of low-mass brown dwarfs up to 30 
Mj, GPI could expect to find at least twice as many 
brown dwarfs as planets. It should quickly become clear 
whether the brown dwarf desert extends to large semi- 
major axes. 

6.2. Model differentiation 

The core accretion model of iMarley et al.l (|2007f ) and 
the hot start model of IBurrows et alf l 20031 ) make dis- 
tinct predictions for the expected properties of young gas 
giants. The differences in formation models are greatest 
for the planets that are young, so a survey targeting stars 
with ages ~ 100 Myr and younger is most effective. If 
gas giants were were formed via only one of the two sce- 
narios, conducting a survey of 100 stars with these ages 
would be sufficient to determine which model described 
the formation of the planets. We performed a simulated 
survey like this twice, once for each model. Figure H] 
shows the expected results of such each survey, with ef- 
fective temperatures of detected planets as a function of 
planet age. It is clear that the two populations of planets 
are very different, both to the eye and from performing 
simple statistical tests comparing the distributions and 
their agreement with each of the models. For this simple 
scenario then, a modest amount of observing time can 
distinguish the two formation models. 

In reality, this dichotomy will not exist. The mod- 
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Fig. 9. — Effective temperatures of detected planets are shown as 
a function of planet age as circles, with representative error bars in 
the lower left corner. Blac k circles indicat e planets with properties 
according to the model of Burrows ct al. (2003), and cyan circles 
correspond to planets from the IMarley et al. (2007) model. The 
lines show expected planet cooling tracks for each model for planet 
masses of 1, 2, 4, 6, 8, and 10 Mj, and the colors f or the cooling 
tracks are black for Burrows ct al. ( 2003) and cyan for Marlcy ct al. 
(2007). More massive planets are at most ages wa r mer t han less 
massive planets in each model, but the Marlcy ct al. (2007) planets 
all converge to roughly the same temperature at young ages. 



els represent two extreme cases, and the post forma- 
tion entropy of any individual planet would likely fall 
somewhere in between what each of the models predicts. 
Looking at a distribution of planets will be even more 
complicated if planets form via more than one mecha- 
nism. It will likely be valuable to segregate distributions 
of detected planets by semi-major axis, host mass, and 
planet metallicity, according to predictions about how 
the planet distribution will depend upon these stellar 
properties for each formation mechanism. Significant dif- 
ferences in these segregated distributions could provide 
useful evidence for a difference in formation history. 

In simulating formation model differentiation, the 
other major simplification was ignoring the existence of 
brown dwarfs, which will also likely complicate efforts 
to determine how planets form. Well constrained stellar 
ages will be vital in discriminating between planet forma- 
tion models. Objects that are brown dwarfs according to 
both the deuterium fusion limit and in the sense of form- 
ing like a star are expected to be significantly warmer and 
more luminous than planet mass objects at the same age. 
Objects with masses above the deuterium burning limit 
that formed via core accretion are potentially more prob- 
lematic if their post formation luminosities are similar to 
planet mass objects that formed through core accretion. 
However, the surface gravities of objects above and be- 
low the deuterium burning limit will be different enough 
to be distinguished, as GPI is expected to achieve a pre- 
cision of better than 0.2 dex in \og(g). 

The HR 8799 system provides an early example of 
examining planet formation models using the effective 
te mperatures and ages of detected planets. As noted 
in iMarois et al.l ([2008a), the members of the HR 8799 
system are not cons istent with th e simp le core accretion 
model presented in IMarley et al.l (|2007l ). After the for- 
mation event, planets formed in their model are insuffi- 
ciently hot to match the temperatures of the HR 8799 
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members, even considering the large uncertainty in the 
age of the planets. Given the dynamical constraints on 
the syste m that place an upper limit on the masses of 
~ lOMj (jFabrvckv fc Murrav-Clavl I2010I ). and despite 
the fact that the thermal history of a planet formed via 
core accretion is und oubtedly more comp licated than the 
model presented by iMarlev et al.l ((2007) , the tempera- 
tures of the HR 8799 planets support the notion that 
planets do not form by core accretion at large semi-major 
axes. The observed abundance of multi-planet systems 
suggests that some of the planets GPI will discover will 
also be in multi-planet systems. Finding more systems 
like HR 8799 will help determine more confidently how 
planets form. 

7. CONCLUSIONS 

We have simulated the performance of the Gemini 
Planet Imager in a variety of hypothetical direct imaging 
surveys of nearby stars, finding how the detection rate 
and properties of detected planets arc affected by differ- 
ent survey choices and assumptions about the exoplanet 
distribution. These simulations rely upon calculations 
and simulations of the noise properties of the system. 
They also rely upon models of the formation and evo- 
lution of gas giant pl anets and their atmo spheres, one 
assuming a hot start ((Burrows et al.l [20031) and one as- 
sumi ng a planet formed by core accretion (jMarlev et al.1 
[20f)l . 

Regardless of formation model, detection rates will be 
highest for observations of young stars. For planets that 
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formed via core accretion, roughly 10% of planets around 
stars younger than 100 Myr may be detected, and for hot 
start planets this detection rate may be as high as 25%. 
One major uncertainty is the frequency of planets at large 
semi-major axes. If planets are very rare beyond 30 AU, a 
GPI survey would be more successful focusing on moder- 
ately young nearby stars, rather than the youngest stars. 
Only by sampling the youngest stars though may GPI be 
expected to place significant constraints on the low-mass 
end of the semi-major axis distribution, as the least mas- 
sive planets cool the most quickly, and sampling young 
stars will also be the most useful approach to better un- 
derstanding how planets form. 

For planets detected with reasonably high signal to 
noise, GPI will be capable of measuring the effective 
temperature and surface gravity. When the age of the 
system is known, these quantities can be compared to 
models of planet formation and evolution. For the ide- 
alized situation where all planets are forme d via either 
the hot start model of iBurrows et all (120031) or the core 
accretion models of IMarlev et al.l 1)20071 ) . a few tens of 
detected planets younger than 100 Myr would be suf- 
ficient to determine the model from which the planets 
were drawn. Real planets will not be so simple, but this 
result suggests that the characterization of many planets 
in this fashion will lead to a better understanding of how 
planets form. 
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